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Abstract. We derive a matrix product formula for symmetric Macdonald 
polynomials. Our results are obtained by constructing polynomial solutions 
of deformed Knizhnik-Zamolodchikov equations, which arise by considering 
representations of the Zamolodchikov-Faddeev and Yang-Baxter algebras in terms 
of t-deformed bosonic operators. These solutions are generalised probabilities for 
particle configurations of the multi-species asymmetric exclusion process, and form 
a basis of the ring of polynomials in n variables whose elements are indexed by 
compositions. For weakly increasing compositions (anti-dominant weights), these basis 
elements coincide with non-symmetric Macdonald polynomials. Our formulas imply 
a natural combinatorial interpretation in terms of solvable lattice models. They also 
imply that normalisations of stationary states of multi-species exclusion processes are 
obtained as Macdonald polynomials at q = 1. 
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Symmetric Macdonald polynomials m 121 are families of multivariable orthogonal 
polynomials indexed by partitions whose coefficients depend rationally on two 
parameters q and t. In the case q = t they degenerate to the more familiar Schur 
functions which encode characters of irreducible representations of the symmetric group. 
One way to define Macdonald polynomials is as joint eigenfunctions of a family of 
commuting operators in the double affine Hecke algebra 011]. They can also be defined 
combinatorially as generating functions nw , or via symmetrisation of non-symmetric 
Macdonald polynomials that are computed from Yang-Baxter graphs 00. 

Another family of commuting operators can be constructed in the Hecke algebra, 
corresponding to commuting transfer matrices in the theory of Yang-Baxter solvable 
lattice models [lO]. Those operators are solutions of the Yang-Baxter algebra generated 
by an i?-matrix arising from a quantum group. Eigenfunctions of transfer matrices 
are generally complicated objects that can be constructed using the so-called Bethe 
ansatz. In cases where the eigenvalue is simple it is sometimes possible to construct 
explicit eigenfunctions of transfer matrices, and a well developed technique to construct 
such eigenfunctions is the matrix product algebra m for the asymmetric exclusion 
process. In an inhomogeneous setting the matrix product algebra is known as the 
Zamolodchikov-Faddeev (ZF) algebra [I2l [13]. The latter algebra is related to the t- 
deformed Knizhnik-Zamolodchikov equation mils]. 

Opdam [16] and Cherednik 00 generalised Macdonalds construction to a non- 
symmetric setting and defined families of non-symmetric Macdonald polynomials which 
are indexed by tuples of integers called compositions. There exists a basis in the 
ring spanned by non-symmetric Macdonald polynomials whose defining equations are 
exchange relations which are equivalent to the t-deformed Zamolodchikov-Faddeev 
algebra. Elements of this basis specialise to probabilities of particle configurations of 
the multi-species asymmetric exclusion process. Solutions to the ZF algebra can be 
obtained from solutions of the Yang-Baxter algebra, but not all solutions of the latter 
give non-trivial solutions of the former. Polynomial solutions of the deformed ZF and 
KZ algebras have received recent attention in the context of the Razumov-Stroganov- 
Cantini-Sportiello theorem and alternating-sign matrices mi nHiiiniEni Ea EH [23] . 

In this paper we put several ingredients together to show that a solution of 
the Yang-Baxter algebra arising for the multi-species asymmetric exclusion process 
and based on t-deformed bosons generates non-trivial polynomial solutions of the 
Zamolodchikov-Faddeev algebra. We show that those solutions lead to matrix product 
formulas for certain basis functions in the ring of polynomials spanned by the non- 
symmetric Macdonald polynomials, and by symmetrisation this implies a matrix product 
formula for symmetric Macdonald polynomials. We furthermore provide a combinatorial 
interpretation of our formulas in terms of two-dimensional lattice configurations of boson 
trajectories. 

For q = 1 the homogeneous limit of our formulas are steady state probabilities of 
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multispecies asymmetric exclusion processes. In the context of the simple asymmetric 
exclusion process the ZF algebra was studied in [2l] and recently reviewed in a 
more general context of integrable exclusion processes in [25]. Our construction here 
is heavily based on recent results for the multispecies asymmetric exclusion process 

[261 EH EH EHl EDI 113]. 


1.1. Deformed bosons 

In the following we shall make extensive use of bosonic operators called t-deformed 
oscillators. These operators appear naturally in the context of solvable lattice models 
and in matrix product formulas for the asymmetric simple exclusion process (ASEP) 
They are dehned by 


ak = tka, 
aa^ = 1 — k, 


a^k = t ^ka\ 
a^a = 1 — t~^k. 


( 1 ) 


A consequence of these relations is the following equation which we shall use frequently, 

aa^ — ta^a = 1 — t. (2) 

A faithful representation of this algebra on the Fock space Span{|m)}5^^o given by 
a|m) = (1 — f’") 2 |m — 1), 


a^|m) = (1 — \m + 1), 


k = diag{l, t, ...}. 
1.2. Hecke algebra action on polynomials 


(3) 


The second algebraic ingredient we need is an action of the Hecke algebra on polynomials. 
We start by denoting by S'„ the symmetric group on n elements, i.e. Sn is the Weyl 
group of type An-i. The symmetric group has a natural action on compositions A G Mq 
given by 

'Si(...) Aj, Aj_|_i,...) (• • •) Aj-|_i, Aj,...). (4) 

The symmetric group also has a natural action on polynomials, given by 


^if {.■ ■ ■ i^ii ^i+1 ) • • •) /'(•••) ^i+1 • • •) • 


(5) 


A t-deformation of ([5]) may be dehned using divided differences leading to the 
Demazure operator 

^±1 ^ ^± 1/2 _ _ s .), ( 6 ) 


^2+1 

which for f —)• 1 reduces to ([5])- The operator Tj satishes the relations of the Hecke 
algebra of type A„_i, 

(T, - =0, (i = 1,..., n - 1) 

rri rri rri rri rri rri 

^ i — 2-^ idil ? 


( 7 ) 
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It is convenient to define the following shifted operator (sometimes referred to as 
Baxterised operator), 

i-l/2 1 — 

Tiiu) = Ti + M = (8) 

which satishes the Yang-Baxter equation, 

Ti{u)Ti+i{u + v)Ti{v) = Ti+i{v)Ti{u + v)Ti+i{u). (9) 

We furthermore dehne the following shift operator compatible with the affine Hecke 
algebra of type A„_i, 


{ujf){Xi, ...,Xn) = f{qXn,Xi, . . . (10) 

ojTi = Ti+iu. ( 11 ) 

The affine Hecke algebra formed by the operators ([H]) and u has an Abelian subalgebra 
generated by the Murphy elements, which are dehned as 

= (12) 

These operators mutually commute and central elements in the Hecke algebra can 
be constructed by taking symmetric combinations of the Murphy elements. Joint 
eigenfunction of the operators Yi are non-symmetric Macdonald polynomials. 

For later use it is convenient to develop some notation. Representations of Sn 
are indexed by partitions, which are those compositions for which Ai > A 2 > ... > A„. 
Compositions are naturally ordered under the dominance order > on fjg, which is dehned 
as 

k 

X>fi if ^^(Aj — Pi) > 0 for A: = 1,..., n. (13) 

i=l 

Let A be a composition and let A"*" be the dominant weight of A in the dominance 
order on compositions, i.e. A"*" is a permutation of A such that Ai > A 2 > ... > A„. 
Let tc_|_ be the smallest word such that A = • A"*". The permutation is obtained 

by labeling each entry from A with a number from 1 to n, from the biggest entry to 
the smallest and from the left to the right. For instance A = (3, 0,4,4, 2) = 

(3, 5,1, 2,4) and so w+ = (3,4,1, 5, 2) and A"*" = (4,4, 3, 2, 0). The anti-dominant weight 
for this example is J = (0, 2, 3,4,4). 

To be able to dehne the spectral vector in the next section we will need the quantity 

P(A) := w+ ■ p, (14) 

where p = |(?7, — 1, n — 3,..., —(n — 1)). For the example above, n = 5 and therefore 
p = (2,1, 0, -1, -2) and p(A) = (0, -2, 2,1, -1). 




Matrix product formula for Macdonald polynomials 
1.3. Deformed Knizhnik-Zamolodchikov equations 


5 


We are interested in finding explicit matrix product formulas for Macdonald 
polynomials. We first derive matrix product formulas for polynomials which are 
linear combinations of non-symmetric Macdonald polynomials. These polynonmials are 
defined by a system of exchange equations closely related to the g-deformed Knizhnik- 
Zamolodchkov {qKZ) equation. We note here that in our notation the parameter that is 
usual called q is replaced by t in order to make connection to the literature on Macdonald 
theory. 

Following Kasatani and Takeyama [20], polynomial solutions to the (reduced) gKZ 
equations, sometimes called exchange relations, can be obtained from eigenfunctions of 
the Yi operators. Let 6 be the anti-dominant weight and let Es be the non-symmetric 
Macdonald polynomial solving the eigenvalue equation 

YiEs = yi{S)Es, (15) 

where yi{S) = Wg now define another set of polynomials, which are linear 

combinations of the Macdonald polynomials Ex, and for which we will be able to find 
compact explicit expression. Define 


fs := Es, 

/...,y,y+i,... := f...x+u\i,... A* > ^*+1- 


(16) 


Note that our notation differs slightly from [20] as fflSj) contains a factor on the 
right hand side, and consequently some other details below are different. 

Then / solves the gKZ equations 


Tif,„x,\ 


i+lv 






1 ^ /...,Ai,Ai+i,... \ ^i+1) 

Ai > Ai+i, 

g^"/Ai,...,A„- 


(17) 

(18) 
(19) 


Writing g = we define the elements of the spectral vector (A) of a composition A as. 


(A)j = Pi{X) uXi. yi{X) = 


( 20 ) 


The non-symmetric Macdonald polynomials are obtained from Es by the action of 
Baxterised operators: 


Esi\ — Ti((A)j_|_i — {X)i)Ex, Xi < Aj+i. (21) 

These polynomials satisfy 

Y,Ex = y^{X)Ex. ( 22 ) 

Below we will derive a matrix product expression for the functions fx, and while 
we will not give an explicit matrix product expression for the polynomials Ex, these 
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can be derived as linear combinations of the fx nsing fl2ip . In fact, the two families of 
polynomials are related via a triangnlar change of basis: 

Ex = Y1 cxM^ f^ = Yl 

/i<A 

for snitable rational coefficients cx^{q,t) and d\^{q,t). 


2. Matrix Product, Yang—Baxter and Zamolodchikov—Faddeev algebras 


2.1. Matrix product Ansatz 

The aim of this section is to obtain a matrix prodnct expression for Es{xi,... ,Xn), 
the non-symmetric Macdonald polynomial indexed by the anti-dominant weight 6. Onr 
approach is to write an Ansatz for the polynomials f\, which generalize E^, and to show 
that this Ansatz obeys the qKZ eqnations (fT7|) - flT^ . The Ansatz is as follows: 


^x+fx{xi, ...,Xn)=TT 


Ax.,{xi) ■ --Ax^MS , 


(24) 


where f2A+ is a normalisation factor to be determined later and Ao(x), Ai(x),..., Ar{x) 
and S are matrices satisfying the exchange relations 

Ai{x)Ai{y) = Ai{y)Ai{x), (25) 

tAj{x)Ai{y) - ^ (^Aj{x)Ai{y) - Aj{y)Ai{x) \ = Ai{x)Aj{y), (26) 

SAi{qx) = q"Ai{x)S, (27) 


for all 0 < z < j < r. It is straightforward to demonstrate that the Ansatz fl2T)) is a 
faithfnl solntion to the qKZ relations flT7)l - (lT^ . The validity of (ITTll is ensnred by (l25ll . 
since fx{xi,..., Xn) is symmetric in Xi, Xi+i when A* = Aj+i. Accordingly, the action of 
(1 — Sj) gives zero in this case, and we hnd that 


'Eifxi,...,Xi,\i+i,...,\ni^ij ■ ■ ■ ! ^n) t ^ f\l,...,Xi+l,\i,...,\ni^ll • • • ! ^n) 


(28) 


by inspection. In a similar vein, when Aj > Aj+i, direct application of fl26|) allows ns to 
conclnde that 


Tif \l,...,Xi,Xi.^.l,...,X„(^^l■! ■ ■ ■ 1 ^n) t /Ai,...,Ai+i,Ai,...,An (^1 ) • • • ) ^n) ■ (^9) 

Finally, to prove 0191) . we observe that 

^A+/A„,Ai,...,A„_i(ga:„,a:i,.. .,Xn-i) = Tr {AxSqXn)MAxi) ■ • • (a^n-i)*^) , 

= Tr (AAi(a^i) • ■ ■ ^x.a-ii.Xn-i)SAxSq^n)) 

= q^"^X+fxi,...,xAxi, (30) 

where we have nsed the cyclicity of the trace and the exchange relation fl27)) to reach 
the hnal eqnality. 
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Solutions to the relations (l25|l - (l26ll can be recovered from the Yang-Baxter algebra 
corresponding to the quantum group or rather a twisted version of it [35] , 

For models based on U£l/ 2 {A^^)., the i?-matrix can be expressed in the form 


r+1 

(x, y) = Y^ 0 ® ® 

i=l ^ l<i<i<r+l 


t - 1 
tx — y 


® ® ( 31 ) 


l< 2 <j<'r+l 


where F'bi) denotes the elementary (r + 1) x (r + 1) matrix with a single non-zero entry 
1 at position {i,j). The intertwining equation, or Yang-Baxter algebra, for such models 
is given by 


y) ■ [Hx) 0 L{y)] = [L{y) 0 L{x)] ■ R{x, y), (32) 


where we have suppressed the superscript (r) and L{x) = LT)(x) is an (r + 1) x (r +1) 
operator-valued matrix. This algebra is well-studied and many solutions for L{x) are 
known. In the next section we will provide a solution for which the elements of L{x) 
are given in terms of f-deformed quantum oscillators. 

The exchange relations fl25l) - fl26l) are equivalent to the Zamolodchikov-Faddeev 
(ZF) algebra [12],[13], 

R{x, y) • [A(x) 0 K{y)] = [K{y) 0 A(x)] , (33) 

where again we suppress (r) and A = AT)(x) is an (r + l)-dimensional operator valued 
column vector given by 

A«(x) = (Ao(a;),...,A,(x))^. (34) 

Equation fl^T]) is instead rewritten as 

SA{qx) = A{x)S, (35) 


where the rank r is again implicit, i.e. A = AT)(x) and S = 

We can construct solutions of fl33l) by rank-reducing the Yang-Baxter algebra fl3^ 
in the following way. Assume a solution of the following modihed RLL relation 


I&\x,y) ■ L{x) ® L{y) = L{y) ® L{x) ■ ^\x^y), 


(36) 


in terms of an (r + 1) xr operator-valued matrix L{x) = and an operator s = 

that satishes 

sL{qx) = L{x)sq~E’‘^^“\ (37) 

which in components just means 


sLij{qx) = q" Mij{x)s. 


( 38 ) 
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(x) = (a:) • (x) • • ■ L^^'> (x), (39) 

sir) ^ ^(r) . ^(r-1) . . . ^(1) 

gives a solution to fl55D and fl5S]) provided that the operator entries of L*^“^ (x) commute 
with those of L^^l{y), for all a ^ b. The usual way to ensure this commutativity is to 
demand that the entries of act on same vector space 14 while act on a different 
vector spaces 14, and indeed we shall adopt this approach in the coming sections. We 
will show that solutions to (l36|l can be obtained from the Yang-Baxter algebra (|3^ by 
trivialising the representation of one of the quantum oscillators, with the consequence 
of reducing the rank of L^'^\x) by one and thus giving rise to L^'^\x). 


3. Low rank examples 

3.1. Rank 1 solution to ZF algebra 

Before presenting the general construction we will first display some explicit examples. 
For convenience we dehne the following functions 


tx-y 

Then for r = 1, we can trivially solve 


t(x - y) 



tx — y ^ 

tx — y 


1 

C--1 

(41) 


tx — y 


for L(^)(x): 


/ 1 

0 

0 

0 \ 

0 

c~ 

&+ 

0 

0 

b- 

c+ 

0 

VO 

0 

0 

1/ 


1 

X 


1 

y 


1 

y 


1 

X 


Hence we see that 


is a rank 1 solution to 
is equal to 


L^^^(x) = A*'^^(x) = 




(42) 


. The corresponding solution to the Yang-Baxter algebra 
1 


L^^\x) = 


xaf X 


(43) 


where the operators a, and k satisfy the f-oscillator relations ([T]). We note that 
trivialising the f-oscillator by sending a\a i—)■ 1 and fc 0, we reduce the rank, and 
thus obtain the solution L(^)(x): 


a 


it 


I-)- 


xa' X 


1 1 
X X 


(44) 
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The rank 2 case is less trivial, giving rise to operator valued solutions for In 

the case r = 2, solving equation (136|1 for we hnd that 


/1 

0 

0 

0 

0 

0 

0 

0 

0 \ 

0 

c~ 

0 

&+ 

0 

0 

0 

0 

0 

0 

0 

c~ 

0 

0 

0 

6+ 

0 

0 

0 

b- 

0 

c+ 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

c~ 

0 


0 

0 

0 

h- 

0 

0 

0 

c+ 

0 

0 

0 

0 

0 

0 

0 

h- 

0 

c+ 

0 

Vo 

0 

0 

0 

0 

0 

0 

0 

1/ 






\( 

1 

a 

\ 


yk 0 

\ya^ y ) 


\ 


7 1 

a\ 

( 1 

a\ 


xk 

0 

$§ yk 

0 


\ xaf 

X j 

V 

y j 



1 a ^ 
xk 0 
xaf X j 


/ 1 

0 

0 

0 \ 

0 

c~ 

6+ 

0 

0 

b- 

c+ 

0 

VO 

0 

0 

1/ 


(45) 


Using (l39i) we thus construct a solution of the ZF algebra in the following way: 


hP‘'^[x) = lP^[x) ■ L^^\x) 


/1 

xk 

a\ 

0 


^ 1 -1- xa ^ 
kx 

yxaf 

X y 

v) 

yxa^ + x^ j 


The associated rank 2 solution to the Yang-Baxter algebra is 


/ 


L^‘^\x) = 


\ 


xa\k2 

xa\ 


ai 

02 

xk2 

0 

xaiol 

X 




/ 


(46) 


(47) 


where and { 02 , 02 ,^ 2 } are two commuting copies of the f-oscillator algebra 

fj47l) . The map aj, Oi 1 and ki ^ 0 reduces the rank of L^‘^\x) by one 



( 1 

1 

02^ 


/ 1 

02^ 

L^‘^\x) !-)■ 

xk2 

xk2 

0 

^ L^‘^\x) = 

xk2 

0 



xa\ 

X y 


\xal 

X y 


(48) 


where the indices of f-oscillators are redundant in the hnal matrix, since we no longer 
need to distinguish between the two copies of the algebra. 

We note that low-rank examples of L-matrices based on deformed oscillators have 
been treated in earlier literature, for example, [211 EH EH EH SO]- 
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We look at an explicit example for rank 2 taking 5 = (0,0,1,1,2,2). In the notation 
of [20] this leads to p{5) = |(—3, —5,1, —1, 5, 3). The nonsymmetric Macdonald 
polynomial that generates solntions to the qKZ eqnation satishes the following eqnations: 


Using the notation 


Y^Es = t-^/^Es 

Y^Es = qt ^Es 


Y 2 ES = t-^/^Es 

Y^Es = qH^'^E, 

(49) 

Y 3 ES = qtEs 

Y^Es = qH^/^Es. 


q = E, 

1 -t- 
^ ^ 1 - f ’ 

(50) 


it can be verihed that the following polynomial solves 


Esi^xi, ...,xe,q,t)= x^x^xp^ + 


[3 + 


{Xi + X2)X3X4X5X6(X5 + Xe) 


+ 


^ 12 ] 


[3 + u][A + u] 


X 1 X 2 X 3 X 4 X 5 XQ. (51) 


We now verify the matrix prodnct form 

^ls+Es{xi,... ,xe-,q = t^,t) = Tr [AQ{xi)AQ{x2)Ai{x3)Ai{xi)A2{x^)A2{xQ)S], (52) 
for this explicit solntion. From (H6|) we see that 

Aq{x) = 1 + xa, 

yli(x) = xk, (53) 

A 2 {x) = xa^ + x^, 

and using (1271) we note that S should satisfy 

qSa — aS = 3, So!' — qaf S = 0. (54) 

Taking the explicit representation ([3]) for the f-oscillators, S has the form 

S = k^ = diag{l, r, ...} = diag{l, q, (55) 


Up to a normalisation, the nonsymmetric Macdonald polynomial E^ is now represented 
in matrix product form by 

Tr [(1 + xio) (1 + X 2 a) x^kxpx^ (a^ + X 5 ) xe (a^ + xe) S~\ 

= X 3 X 4 X 5 X 6 Tr [(xsXgA;^ + (xi + X 2 )(x 5 + xe)ak‘^a^ + xiX 2 a^A:^(a^)^) S~\ , (56) 

where other terms involving unequal powers of a and of have zero trace. Normalising 
with = 14221100 = Tr(A:^S') we hnally get 


Es{xi, ...,XG;q = r,t) 


X 3 X 4 XPI + X3X4X5Xq{Xi + X2)(X5 + Xe)^ 


.Ti aa^k^S 


+ XiX2X3X4X5X6t' 


Tr k^S 
4Tra2(at)2pF 

TrFF ’ 


(57) 


which can be shown to equal fl5T|) . We give the details of calculating the traces in 
Section [HI 
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The ZF algebra for the case x = ?/ = 1 is known as the matrix product algebra for 
the multi-species asymmetric exclusion process. A general rank solution for this case, 
hrst for t = oo and then in terms of f-oscillators, was recently obtained by several 
authors in a sequence of works [26l [23 [2H1 EHl [30]. A generalisation of these results 
that includes a spectral parameter which can be chosen inhomogeneously was found 
earlier in [31] and independently in the case of super-algebras in [32] . In Section 14.11 
we shall illuminate the very natural combinatorial structure of this solution. We note 
that a different inhomogeneous generalisation of the multi-secpies ASEP, with species- 
dependent hopping parameters, was studied in [H] HS] US] Hi] . 


dr). 


(x) whose entries 

are given by 


^ nm=j+l ^rn, 

i = j 


Xttjttl nm=i+l ^rn, 

i > j 

(58) 

0 , 

i < j 



for all 1 < i, j < r, and 


r 

= Oj, I < j <r, L^ioix) = xa\ km, I < i < r, 


4 ? = 1 . ( 59 ) 


where {ai,al,ki}, 1 < i < r are r commuting copies of the t-oscillator 
Then this L matrix satisfies the intertwining eguation fl32l) . 


algebra O^. 


An example of of L for r = 3 is given in flTOj) . A proof of Theorem [T] can be obtained 
by a long brute force check of the intertwining equations fl32l) . by distinguishing many 
different cases. Here we present a more compact and elegant proof. First we start with 
an easy general property of the solutions of the intertwining equations (l32ll . 


Lemma 1. Let L{x) be a {r + 1) x {r + 1) matrix of operators depending on a spectral 
parameter x, and {no,...,^^} and {uq, ... ,Ur} two constant vectors of commuting 
operators that also commute with the entries of L{x). Define the matrix 


C{x) = diag{no,..., n J ■ L(x) • diagl^o, 


(i) If L{x) is a solution of the intertwining eguations then C{x) solves 

(a) Conversely, if C{x) is a solution of the intertwining eguations and the operators 
Vi,Uj are non identically zero, then L[x) solves / f^) . 

Next we want to describe how to get a solution of the intertwining equations fl32l) of 
ra nk r from a solution of rank r — 1. For 0 < m < r — 1 introduce an insertion operator 
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Mij 

0 < z < m — 1, 

0 <j<r— 

1 


i = m, 

0 < j < r, 


Mi—ij 

m + 1 < i < r, 

0 < j < r — 

1 

0 

0 < i < m 

j = r. 



Here is an example for r = 3, 


9 = 


^ 9oo 901 9 o 2 ^ 
9io 9n 9 i2 

y 920 921 922 j 




O 

O 

9oi 

902 

0 ^ 

0 

0 

0 

n 

9io 

911 

912 

0 

\ 920 

921 

922 

0 ^ 


(60) 


(61) 


Lemma 2. Let L{x) be a solution of the intertwining equations of rank r — 1 such 
that the first m rows are independent of x, while the last r — m are linear in x. The 
matrix Xm,K {.L{x)) is a solution of the intertwining equations / El) of rank r if and only 
if the operator n satisfies the following commutations 


KLij{x) = Lij{x)K Vj, 0 < z < m — 1 
KLij{x) = tLij{x)n Vj, m < i < r — 1 


(62) 


Proof. The proof of this lemma is a straightforward check of (l32|) . 


□ 


We have already seen that the matrix 

L{a,a\x)= \ ^ I (63) 

\ xa' X I 

satishes the intertwining eqnations fl5^ of rank 1, and we set 

L^^\x) = L{ai,a\,x). (64) 

Now we are going to constrnct L^'^\x) starting from in a reenrsive way. The 

hrst ingredient we need is the anxiliary matrix 


= \Po (diag{l, Qo, • • •, Qo} ■ 


where Pq, Qq commnte with the entries of ^)(a:) and satisfy 


(65) 


PqQo — tQoPo- (66) 

Assnming that satishes the intertwining eqnations and applying Lemma [Hand 

Lemma [2] it follows that £)’''’^)(x) satishes the intertwining eqnations fl32l) . 

The second ingredient we need is the matrix jC[j’''\x) dehned by 

=X^_^^P^kr-- ■A,Pikr (L (o^, 4, x) •diag{l,Qi}) 


(67) 








Matrix product formula for Macdonald polynomials 

where here Pi, Qi commute with a^, aj, kr and t-commute among themselves 


13 


PiQi — tQiPi- 


By applying Lemma [Hand repeatedly Lemma [2] we have that also satisfies the 

intertwining equations fl3^ . 


Then we put the two ingredients together and define 


£h+h(cc) ;= C^+P{x) ■ C^+P{x) 


( 68 ) 


which obviously satisfies the intertwining equations fl32l) if L^'^\x) does. Comparing 
£(r+i)(^) we see that 

£h+h(a:) = diag{l, QoA, QoA, • • •, QoA, 1} ■ {x) ■ diag{l,..., 1 , QiPq}, (69) 

therefore Theorem [1] follows from the second point of Lemma [Hand induction on r. 


4-1- Graphical interpretation of L-matrix 

While the L-matrix as given by equations f[58|) and f[5^ solves the intertwining equation 
(j32|l . the precise form of the entries is quite mysterious without any further explanation. 
The aim of this section is to show that all entries can be deduced from some governing 
combinatorial rules. To that end, we denote the entries of the L-matrix graphically by 
tiles. The left and right edges of the tiles are either unoccupied (corresponding with 
index 0), or occupied by a boson of colour i (corresponding with index i). For example, 
for rank 3, the entries of the L-matrix are encoded as follows: 



where the indexing conventions are o = 1, o = 2, • = 3, with an empty edge representing 
the index 0. 

This information alone is sufficient to label the entries of the L-matrix 
unambiguously. But for the purpose of motivating the form of the entries it is useful to 
consider, further to this, lattice paths which are generated by the bosons. The horizontal 
edges of tiles can be occupied by arbitrary numbers of bosons from each of the families, 
and we adopt the convention that the bosons are ordered from darkest to lightest, 
reading from left to right. In passing from the bottom horizontal edge to the top one, 
the occupation number of any species of boson can go up/down by 1 (representing the 
action of a creation/annihilation operator), or remain the same (representing the action 
of some diagonal operator). We keep track of these transitions between states by simply 
“connecting the dots”. For example, 
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(71) 


represents the action of the operator a^Oi (creation of a type 2 boson, annihilation of a 
type 1, with time flow up the page). These graphical conventions suffice to explain all 
creation/annhilation operators appearing in the L-matrix (JTOl). 

In order to specify the zero entries of flTOjl . as well as the inclusion of the ki operators, 
we need two rules for the lattice paths thus constructed: 

1. It is forbidden to create a boson of type i and simultaneously annihilate one of type 
j, when i < j. 

2. We obtain a factor of t every time a type j line horizontally crosses a type i line, 
where i > j. 

Taking into account these rules, we can now reproduce the exact form of all entries 
of the L-matrix. The zero entries correspond with tiles forbidden by rule 1. Rule 2 
explains which entries are dressed by ki operators. For example, (na corresponds to 
L^l = k^alai, with k^ producing a factor of t^^, where m 3 is the number of type 3 bosons 
present. We easily deduce this from rule 2, since we will have exactly m 3 horizontal 
crossings of type 3 lines by the left-turning type 2 line. A further example: 



(3) t 

corresponds with L) q = k 3 k 2 a[, since the left-turning type 1 boson horizontally crosses 
all type 2 and 3 lines present, and we must keep track of these crossings by the inclusion 
of k3k2. 


4 . 2 . Trivialising a representation 

To obtain a solution to the rank-reduced intertwining equation fl3^ we choose the 
representation of the hrst t-oscillator algebra to be trivial: 

ai = a| = 1, ki = 0. 

With this choice of representation, it is easy to see that the hrst two columns of the 
L-matrix flS5]l . become equal, leading to an immediate solution of fl5B]) (by simply 
omitting one of the redundant columns from the L-matrix). For the example flTOj) , after 
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trivializing the first t-oscillator algebra, we obtain 

02 O3 

xk^k2 0 0 

xk'ia^ xks 0 
xal xala 2 x y 

Up to some elementary transformations, this is exactly the rank 3 solution obtained in 
[29] (see equation (47) therein), although here we write our operators with subscripts to 
distinguish commuting copies of the f-oscillator algebra, rather than the tensor product 
notation employed in [29] , 



L^^\x) = 


/□ □ 
ODD 
D Q Q 
□ DO 


4-3. Solution of ZF algebra 


Having obtained solutions of fl36|) for all r, we construct solutions of the ZF algebra 
via the prescription fl39l) . This is conceptually straightforward, although it introduces 
a slight notational complexity: we must now distinguish not only between different 
families of f-bosons (which we have done so far by using subscripts), but also between 
operators which act at different lattice sites (which we now do by placing superscripts 
on our operators). For example the rank 3 solution of the ZF algebra is given by 


A^^^(a;) = 


/ 

1 

a 2 

as 

\ 

(3) 

/ 

1 


2(^3/C2 

0 

0 




L 


xk^a\ 

xkz 

0 




xk 

V 

xa\ 

xa\a 2 

X 



V 

xa 


a2 \ 
0 


( 2 ) 


X 


( 1 ) 


/ 


/ Ao(x)\ 

Ai{x 

Mx. 

\Mx)J 


(73) 


where the superscript placed on a matrix indicates that all operators within that matrix 
acquire that superscript. In terms of the graphical conventions that we have introduced, 
the components Ai[x) of the rank r solution are given by rows of tiles of length 

r. The left boundary of the row is occupied by a particle of colour f, for 1 < z < r, or 
unoccupied in the case i = 0. The right boundary is always unoccupied. For example, 
formulating fl73|l in terms of tiles, we obtain 

(3) 

/ 


/ 


A‘^^\x) = 


□ 




O D 
O D 
□ D 


□ \ 
D 
Q 
O/ 



f Ao(a:)\ 
Ai{x 
A 2 


Uc 


( 74 ) 


\Mx)/ 


From this it is easy to extract individual components, for example: 


A 


+ 


+ 


2 X = 




(3) (2) (1) 

rk^3) f(3) 

JU ri/Q Lho 


(3) (2) 

™2l( 3) jL(2) 
Jb A/O njry 


( 1 ) 


X 


(3) (2) (1) 
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where the top and bottom edges of the row can be occnpied by arbitrary nnmbers of 
“backgronnd” particles, which we do not show here, since they play no role in specifying 
the operators at each site. 

In the discnssion in Sections 14.II and 14.21 we made no distinction between the colonr 
of a particle, and the family of bosons that it represents. Indeed, a particle of colonr 
i corresponded with the t-oscillator algebra {aj,a|,/cj}. However (with the conventions 
that we have adopted) in passing to the fnll solntion of the ZF algebra, which is valned 
on ( 8 ) • • • ( 8 ) the notions of colour and family no longer coincide across all 
For example in fTM)) . a particle of colonr • corresponds with the algebra { 03 , 03 ,^ 3 } 
when in colnmn 3, with { 02 , al, ^ 2 } when in colnmn 2, and with the trivialized algebra 
{oi, a|, A:i} when in colnmn 1. The general rnle is that a particle of colonr i in colnmn 
r —j corresponds with the f-oscillator family i—j, and this shift shonld be kept in mind 
in the snbseqnent sections. 

f-f - Representation of twist operator 

Having constrncted solntions of (l36|) for all valnes of r, we now seek an explicit operator 
s which satishes fl38l) . This is very easily achieved by considering the form of the entries 
of L{x). We let s be factorized over the non-trivial copies of the t-oscillator algebra, 
with the following commntation relations: 

s = Sr---S2, OiSi = a|sj = kiSi = Siki, (75) 

where as nsnal all operators carry the snperscript (r), which we have snppressed for 
visnal clarity. 

Proposition 1. Equation (|38il is satisfied with s as defined by (|75il . 

We obtain a valid solntion of fl75l) by letting 

Si = = diag{l, ■ ■ .}i (76) 

which is the representation that we nse in all snbseqnent sections. 


5. Transition formulas 


The aim of this section is to write a recnrsion relation for the fnnction f\{xi,... ,x„) 
for any composition A, nsing the fact that both Ah)(x) and are recnrsively dehned. 
Indeed, we can rewrite eqnations fl39|) and fHOj) as 

(77) 

Sir) ^ gir)sir-l) (73) 

where the action in the vector space is explicitly factorized ont. This snggests the 
dehnition of a transition matrix: 




Tr 


(^1) ■ ■ ■ {^n)s 
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where all operators carry a subscript (r) which we have suppressed, and A, n are 
compositions whose parts lie in [0,1,..., r] and [0,1,..., r — 1], respectively. We will 
use this transition matrix to derive the recursion relation that is of interest to us, but 
before we do that, we introduce a simple operation on compositions. Let A* be the 
composition obtained by subtracting 1 from all non-zero parts of A: 

A* = max(Ai — 1,0). 


Theorem 2. Let X be a composition whose parts satisfy 0 < Aj < r, for all 1 < i < n. 
The following recursion relation holds: 

pnj{\) 

where mi{X) is the number of parts of size i in X and the sum is taken over all 
compositions p that are a permutation of X*. 

Proof. Rewrite fl2^ as 

kl\+ /a(xi, ..., Xn) = Tr [^Ai(a^i) • • • , 

for some proportionality factor r2A+ = ^2^+ that remains to be determined (and in 
particular needs to be shown different from zero). We remark that the factor 17^+ is 
the same for all partitions A with the common re-ordering A"*". Combining flHOD with 
equations (1771) and (1781) . we obtain 


Tx,f,{Xi,...,Xn)f^{Xi,...,Xn), (79) 


r-l / 

fx{xi,...,xn )=n 

2=1 \ j<i 




5^ Ti' hi’-y, (,T.)... ij-y, IV [<-‘>(*1) ■.. . (so) 


which can be rewritten as 


kl^^lfxixi, ...,Xn) = ^Tx,^{xu . . .,Xn)kl^^+ . . . ,X„). (81) 

At this stage the sum in (IHT]) is taken over all compositions //, with parts in 
[0, 1,... , r — 1]. In order to reduce the size of this sum, we need the following result. 

Proposition 2. For any two compositions X and p., if ^ (-^*)''' then 

rA,^(xi,... ,x„) = 0 . 

Proof. Given a composition A, let mj(A) be the number of its parts of size i. We wish 
to show that if mj(A) 7 ^ mi-i{p) for some i > 2, then Ta,^(xi, ..., x„) = 0 . 

For each i > 2, let ffi{X,p) be the number of pairs {Xk, Pk) = [ip — !)• Invoking 
the explicit form of L{x), we see that the number of a| operators appearing in 
Tx,k.{^i, ■ ■ ■ ,Xn) is equal to mi{X) — ifi{X,p). Conversely, the number of a* operators 
appearing in this trace is given by mj_i(/z) — ifi{X,p). In order for the trace to be 
non-vanishing, these two numbers must be equal for all i > 2. □ 
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An immediate consequence of this proposition is that we can rewrite 0811) as 

n^^lfx{Xi,...,Xn) ^ Tx,f,{Xi,...,Xn)ff,{Xi,...,Xn) (82) 

MG 5 „-(A*)+ 


To complete the proof of 0791) . it remains only to demonstrate the following result. 
Proposition 3. The proportionality factors appearing in 0821) satisfy 


r—1 



1 

qifK-K+i ’ 


(83) 


where is the conjugate of A+. 


Proof. Without losing generality, we ease notation by assuming that A is a partition. We 
compare coefficients of the monomial x^ := x\^ ■ ■ -xfj" on both sides of 08^ . Since the 
matrix elements Ta,^(xi, ... ,x„) are at most degree 1 in each Xj, the only contribution 
to x^ in the right hand side comes from p = X*. Therefore we can write 




nx’x^ = n 


A. Tx,x*(,^l^ ■ ■ ■,Xn)x^* = ^^Ta,A‘(1, ■ ■ ■ , l)x^ 


(84) 


or more simply, ^^Ta,a»(1 , • • •, !)• A simple, explicit computation gives 

... 1 ) = niv = n 

i=l i=l ^ 

□ 

Thanks to this proposition, flH^ implies fl79|) provided that A^a^T^^ is non-vanishing. 
This can be deduced inductively on r, using flS^ with = 1. Indeed, we find that 



nni 


1 


n 


1 


1 




□ 


A simple case of fl79l) is when all parts of the composition A are non-zero. A* > 0 
for all 1 < z < n. In this case the summation on the right hand side is trivial, and we 
recover 


/a(xi, ..., Xn ) 


JJa;* /a*(xi, .. ., Xn ). 


vi=l 
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6.1. Lattice formulation of matrix product 

In view of the graphical representation of the elements of the L-matrices, it is possible 
to interpret the matrix product formula fl2^ entirely in terms of lattice paths. This 
leads to a combinatorial rule for i Xn), for any composition A. 

We consider a lattice formed by n rows of tiles, with associated spectral parameters 
Xi,... ,Xn, where the index of the spectral parameters increases as we go from the top 
row to the bottom. Each row is r tiles wide, where r is the size of the largest part 
in A. For each A* > 0, there is a boson of colour i incident on the left boundary 
of the row. However, the rows i for which Aj = 0 have no particle at the left 
boundary. The right boundary of the lattice is completely unoccupied by particles. The 
lattice thus constructed reproduces all terms in the matrix product fl2T|) . since each row 
corresponds with one of the components Ai{x) of A^^\x). For example, for r = 3 and 
A = (0, 2,3,1, 0, 2), the matrix product can be represented in the following way: 

Xi 

X2 

X3 

Ti{AQ{xi)A2{x2)A3{x3)Ai{xi)Ao{x^)A2{xQ)S) = X4 

X5 

Xq 


Some care needs to be taken in regard to the boundary conditions on the top and 
bottom edges of the lattice. Since we are taking a trace in fl2^ . the correct way to 
view the lattice is with the top and bottom edges identified, so that lattice paths are 
permitted to “wrap around” the cylinder thus formed. This ensures that in each column 
of the lattice, there is conservation between the number of particles entering/leaving that 
column (in other words, all creation/annihilation operators come in pairs, ensuring a 
non-vanishing trace). The sole exception to this rule involves the boson species whose 
representation has been trivialized. Recall that in column j of the lattice, bosons from 
any of the families {1 ,..., j} are allowed, but we trivialize the representation of the 
first family. Hence at the top of each column we are permitted to eject type 1 particles 
without causing the trace to vanish. In the example above, we eject 1 trivial boson from 
column 3, 2 trivial bosons from column 2, and 1 trivial boson from column 1, since these 
are the corresponding multiplicities in the composition A. 

Having set up the boundary conditions of the lattice in this way, we obtain all 
possible terms in the matrix product fl2^ by summing over all lattice paths which evolve 



( 3 ) ( 2 ) ( 1 ) 
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from the left edge to the top boundary (with round-the-cylinder contributions being 
allowed). Note that each conhguration gives rise to an explicit product of operators, 
whose trace must then be calculated, rather than local Boltzmann weights, as would 
normally be the case in combinatorial calculations of this nature. 

Let us remark hnally that equation fITOl) has a particularly natural meaning in the 
above combinatorial framework. If we represent fx{xi ,..., Xn) in terms of an nxr lattice 
as described, then fl7^ can be interpreted as summing over all possible conhgurations of 
column r of the lattice, with the remaining r — 1 columns giving rise to ffi{xi ,..., Xn)- 
The elements Tx^^{xi,... ,Xn) of the transition matrix can thus be viewed as partition 
functions of a single column. We illustrate the calculation of Tx^^{xi,... ,Xn) on an 
explicit example in the Appendix. 


6.2. Rank 2 example 


We give an example of the lattice-path calculation for fs{xi,... ,Xn), where 6 = 
(0, 0,1,1, 2, 2) is an anti-dominant weight. This is of particular interest, since in this 
case fs = Es{xi, ..., Xn] q, t). The possible lattice conhgurations for this example: 


Xi 

X2 

X3 C 

Xi C 

Tr- i 









) 




) 







nr.n i 










( 2 ) ( 1 ) 


Xi 

X 2 
X3 

X4^ 

X5 

Xq 

(2) (1) 









— H 




d 





o 















( 2 ) ( 1 ) 


—a 

TTl 



o 


o 


r 

-K—1 

j 



(2) (1) 



where we highlight in green the tiles which give rise to an x weight (these are the tiles 
which have an occupied left edge). Note that in column 2 of the lattice it is possible 
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for the type 2 bosons to go around the cylinder, but this does not lead to any particle 
conhgurations not already shown here. From these lattice conhgurations we can read 
off the product of operators which we need to trace over. We obtain 

ils+Es{xi,... ,XQ;q,t) = x^x^x^xl Tr[/c^S'] + X2X3X4X5XI Tr[a/c^a^S'] + XiX'j,X/iX^xl Ti[ak‘^a) S] 
+X2X'iXixlxQ Tr[aA;^a^S'] + xix^x^xIxq Ti[ak‘^a)S] + XiX2X‘iXiX^XQ Tr[aafc^a^a^S'], 

in agreement with our earlier calculation fl57|) . 


7. Symmetric Macdonald polynomials 


We relate our results to symmetric Macdonald polynomials. The following result is 
already mentioned in 


Lemma 3. Let X be a dominant composition, i.e. a partition. Then the sum 

Vxixi,...,Xn;q,t) = ^fi,{xi,...,xn;q,t), (85) 


where the sum runs through all permutations fi 0 /A+, is symmetric. 

Proof. We need to show that TfPx = t^^'^Vx for alH = 1 ,..., n — 1. From flTS]) and 
we hnd for Aj < Aj+i that 

= ( (1 + (*'''" - 


Combining this with fll7p and fllSp we thus hnd 

fi /i: fii<fii+i /i: fi: 

= XZ 'Ll */c = *5Z/c' 


( 86 ) 


(87) 


/i: 


At: 


□ 


The Macdonald polynomial Px is the unique symmetric polynomial (up to 
normalisation) which can be obtained by taking linear combinations of the non- 
symmetric Macdonald polynomials where /r is a permutation of A. By (l2^ . each 
can be written as a linear combination of non-symmetric Macdonald polynomials 
where is a permutation of // and z/ < /r. It follows that Vx is a linear combination of 
such that /i is a permutation of A, and since Vx is symmetric, Vx and Px must be 
equal up to normalisation, due to uniqueness. By this argument and Lemma [3l we have 
the following theorem. 
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Theorem 3. For A C r"' 

n 

QxP\{xi,... ,Xn]q,t) = ^ (88) 

2=1 

where 12^ = S = and the sum is over all permutations p of X. 

A corollary of Theorem [3] is that the generating fnnction for symmetric Macdonald 
polynomials for partitions A C r” can be written as a single matrix prodnct. Define first 

r 

A{yo,...,yr-,x) = ^yiAi{x). (89) 

1=0 

Notice that [A{yo,... ,yr;x), A{yo, ■ ■ ■ ,yr', = 0 by a similar argnment as in Lemma El 

Moreover we have that 

n r 

TT[SYlA{yo,...,yr-,Xi)] = ^ Da JJ ..., g, f), (90) 

i=l ACr" i=0 

where again all snperscripts (r) are snppressed and mi{X) is the multiplicity of part i in 

A. 

Special values of q include the Hall-Littlewood polynomials {q = 0), which have 
been expressed as Bethe wave-functions in the t-boson model [15]. The structure of the 
expression in |15] has some features in common with our matrix product formula, but is 
distinguished by some important differences. In particular, higher rank L-matrices are 
not needed in the Hall-Littlewood formula, and the construction does not involve a trace 
operation. Other special values of q include Schur functions {q = t), Jack polynomials 
{q = t°‘,t ^ 1) and g = 1 results in the normalisation of the inhomogeneous multi-species 
asymmetric exclusion process with periodic boundary conditions. The homogeneous 
limit of the latter is recovered by sending Xj ^ 1. 

8 . Calculation of traces 

So far we have explained how to calculate the matrix product fl2T)) in terms of traces of 
f-oscillators, without explicitly evaluating these. We now turn to the evaluation of the 
traces, using the representation ([3]) of the t-oscillator algebras. Let us start with the 
examples encountered in equation (1571) : 

oo 

TrikP] = = — = 7—TFT- 

Tr[/i;P] Tr[A:P] [p -|- 1] ’ 

TT[aaa^a^kP] _ _ [2] 

Tr[kP] Tr[/cP] \p + l][p + 2]' 


( 91 ) 
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In theory, the trace of any product of f-oscillators can be evaluated by simply summing 
up geometric series. On the other hand, as the examples become more complicated, the 
answer is no longer neatly factorized. In order to evaluate the most general case, let 
be the set of all Dyck paths of length 2i. It is easy to verify that all types of traces 
arising in this work have the generic form Tr[D£A;^] where p is an arbitrary exponent, 
and represents a Dyck path of length 2i generated by the oscillators (a means an 
up step, means a down step). Indeed, the trace of any string of t-oscillators can be 
brought into this form (up to overall factors of t) using the cyclicity of the trace and the 
commutation relations ([T]). Let us denote D^ by a series of open/closed parentheses, for 
example: 

Tr[()(())A:^] := Tr[aa^aaa^a^A:^]. 


We consider the following map: 

M: ^ 

D^ iH- (mi,.. .,mi) 


(92) 


where m^ is the number of parenthetic pairs (...) in D^ which are surrounded by z — 1 
other parenthetic pairs. For example in the case (()(())) we have mi = 1, m 2 = 2 and 
m 3 = 1. Using this map we deduce the following lemma. 

Lemma 4. 

00 

TTpekP] = ... (1 - (93) 

n=0 

where (mi,..., m^) = M. (D^). 

It is now useful to introduce the following operators on functions: 

= fiz) - f{tz), = 2:0 6 to-yo 5t , m > 1, = 2;. (94) 

m 

The operator A^™^^ acts multiplicatively on monomials and it is easy to see that 

= (95) 


It therefore follows that 

CXD 

... (1 - 

00 

^(rne) _ _ _ ^("^ 1 ) 


n=0 




n=0 


J z=l 


^i^e) . . . ^(mi) _^ 


1 — XZ 


(96) 


Z = 1 


If we dehne the function 




(mi)_ 

1 — XZ 


(97) 


J Z=1 
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Tr[I)ikP] = where {mi,... ,mi) = A4(D^). 

Hence we are able to evaluate all traces that we encounter in terms of the function (j97ll . 
whose dehnition is relatively elementary. 

9. Conclusion 

The main result of this paper is a matrix product formula for Macdonald polynomials 
in terms of deformed bosonic operators. This formula implies a new explicit way to 
efficiently compute these polynomials and also provides a combinatorial interpretation. 
Our result hrmly connects the polynomial representation theory of the affine Hecke 
algebra with the theory of solvable lattice models using tools such as the Yang-Baxter 
equation, /^-matrices, the Zamolodchikov-Faddeev algebra and the deformed Knizhnik- 
Zamolodchikov equation. 

The results discussed in this paper have a direct generalisation to the 
inhomogeneous multi-species asymmetric exclusion process with boundaries |1H1 SHI SEl 
Wf\ l50] resulting in a matrix product formula for Koornwinder polynomials |5T]. A 
connection between Koornwinder polynomials and the quantum XXZ chain, which is 
closely related to the exclusion process by a similarity transformation, was made in |52j . 

We note that a different generalisation of the multi-species asymmetric exclusion 
process, using inhomogeneous hopping parameters, was discussed in [H] HSj 021 01] . It 
would be interesting to clarify the matrix product structure |13| in this case, as well as 
the relation to integrability [53] . 

Acknowledgment 

LC has been supported by the CNRS through a Chaire d’Excellence. JdG and MW 
are generously supported by the Australian Research Council (ARC) and the ARC 
Centre of Excellence for Mathematical and Statistical Frontiers (ACEMS). We thank 
Kayed A1 Qasemi, Eric Ragoucy, Sergey Sergeev, Ole Warnaar and Paul Zinn-Justin for 
discussion. 

Appendix A. A further combinatorial example 

We illustrate the combinatorial meaning of on a small rank 3 example, namely 
/a(xi, X 2 , X 3 , X 4 ) with A = (3, 1 , 0 , 2 ). The function /( 3 ,i,o, 2 ) can be expanded into 
functions where is a permuation of (2,1,0,0). So in principle there are 12 terms 
in this expansion but since the entries ( 1 , 1 ), ( 1 , 2 ), and ( 2 , 2 ) of vanish, we have 
only to consider the permutations (tti, 7 r 2 , 713 , 714 ) of ( 2 , 1 , 0 , 0 ) such that 772 7 ^ 1,2 and 
774 7 ^ 2. Hence we obtain the following four terms in the expansion: 

/a = P'a,(2,0,0,1)/(2,0,0,1) + 7a,(0,0,2,1)/(0,0,2,1) + P'a,(2,0,1,0)/(2,0,1,0) + (1,0,2,0)/(1,0,2,0) • (A-1) 
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These terms are represented by the following four pictures: 
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Replacing each tile in column 3 with the operator(s) it represents, this becomes 

/(3,1,0,2) = (1 - qt){l - qH^)xiX 2 Xi (TT[k3k2hs]f(2,o,o,i) + Tr[4A:3A:2a3fc3s]/(0, 0,2,1) 
+ Tr[/i;3A:2a2fc34'S]/(2,0,1,0) + Tr[4a2/c3A:2a3fc34'S]/(i,o,2,o)) • 

The traces can then be factorized across commuting families of operators: 


/(3,1,0,2) = (1 - qt){l - qH^)xiX 2 Xi (Tr[4s3]Tr[A:2S2]/{2,o,o,i) 

+ q^t Tr[a34/c3S3]Tr[/c2S2]/(o,o,2,i) 

+ Tr[4s3]Tr[a24A:2S2]/(2,o,i,o) + gV Tr[a34A:3S3]Tr[a24/c2S2]/(i,o,2,o)) , 
and using the explicit form fl75|) of the twist operators, this becomes 
/(3,1,0,2) = (1 - ?^)(1 - qH‘^)xiX2Xi (Tr[/c^+2“]Tr[/i;^+“]/(2,0,0,1) 

+ q^t Tr[aa'f/c^+^“]Tr[A;^+“]/(o, 0 , 2 , 1 ) 

+ Tr[/c^+^“]Tr[aa'''/i;^+“]/(2,o,i,o) + Tr[aa^/i;^+^“]Tr[aa'f/c^+“]/(i,o,2,o)) • 

The traces themselves are easily calculated using the results in Section |8l Thus we 
obtain the matrix elements rA,(2,o,o,i), rA,(o,o,2,i), Ta,( 2 ,o,i,o) and Ta, (i,0,2,0)- 
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